Abundant incongruence in a clade endemic to a biodiversity hotspot: Phylogenetics of the scrub mint clade (Lamiaceae).

The Scrub Mint clade(Lamiaceae) provides a unique system for investigating the evolutionary processes driving diversification in the North American Coastal Plain from both a systematic and biogeographic context. The clade comprisesDicerandra, Conradina, Piloblephis, Stachydeoma, and four species of the broadly defined genus Clinopodium(Mentheae; Lamiaceae), almost all of which are endemic to the North American Eastern Coastal Plain. Most species of this clade are threatened or endangered and restricted to sandhill or a mosaic of scrub habitats. We analyzed relationships in this clade to understand the evolution of the group and identify evolutionary mechanisms acting on the clade, with important implications for conservation. We used a target-capture method to sequence and analyze 238 nuclear loci across all species of scrub mints, reconstructed the phylogeny, and calculated gene tree concordance, gene tree estimation error, and reticulation indices for every node in the tree using ML methods. Phylogenetic networks were used to determine reticulation events. Our nuclear phylogenetic estimates were consistent with previous results, while greatly increasing the robustness of taxon sampling. The phylogeny resolved the full relationship between Dicerandra and Conradina and the less-studied members of the clade (Piloblephis, Stachydeoma, Clinopodium spp.). We found hotspots of gene tree discordance and reticulation throughout the tree, especially in perennial Dicerandra. Several instances of reticulation events were uncovered between annual and perennial Dicerandra, and within the Conradina + allies clade. Incomplete lineage sorting also likely contributed to phylogenetic discordance. These results clarify phylogenetic relationships in the clade and provide insight on important evolutionary drivers in the clade, such as hybridization. General relationships in the group were confirmed, while the large amount of gene tree discordance is likely due to reticulation across the phylogeny.

[Identification and analysis of terpene synthase (TPS) gene family in Schizonepeta tenuifolia].

Terpenoids are important secondary metabolites of plants that possess both pharmacological activity and economic value. Terpene synthases(TPSs) are key enzymes in the synthesis process of terpenoids. In order to investigate the TPS gene family members and their potential functions in Schizonepeta tenuifolia, this study conducted a systematic analysis of the TPS gene family of S. tenuifolia based on the whole genome data of S. tenuifolia using bioinformatics methods. The results revealed 57 StTPS members identified from the genome database of S. tenuifolia. The StTPS family members encoded 285-819 amino acids, with protein molecular weights ranging from 32.75 to 94.11 kDa, all of which were hydrophilic proteins. The StTPS family members were mainly distributed in the cytoplasm and chloroplasts, exhibiting a random and uneven physical localization pattern. Phylogenetic analysis showed that the StTPS genes family were divided into six subgroups, mainly belonging to the TPS-a and TPS-b subfamilies. Promoter analysis predicted that the TPS gene family members could respond to various stressors such as light, abscisic acid, and methyl jasmonate(MeJA). Transcriptome data analysis revealed that most of the TPS genes were expressed in the roots of S. tenuifolia, and qRT-PCR analysis was conducted on genes with high expression in leaves and low expression in roots. Through the analysis of the TPS gene family of S. tenuifolia, this study identified StTPS5, StTPS18, StTPS32, and StTPS45 as potential genes involved in sesquiterpene synthesis of S. tenuifolia. StTPS45 was cloned for the construction of an prokaryotic expression vector, providing a reference for further investigation of the function and role of the TPS gene family in sesquiterpene synthesis.

Full-length transcriptome sequencing provides new insights into the complexity of flavonoid biosynthesis in Glechoma longituba.

Glechoma longituba has been frequently used in treating urolithiasis and cholelithiasis due to the presence of flavonoids, which are its major bioactive constituents. However, research on the molecular background of flavonoid biosynthesis in G. longituba is limited. In this study, we used single-molecule real-time combined with next-generation sequencing technologies to construct the complete transcriptome of G. longituba. We identified 404,648 non-redundant transcripts, including 249,697 coding sequences, 197,811 simple sequence repeats, 174,846 long noncoding RNA, and 176,554 coding RNA. Moreover, we functionally annotated 346,218 isoforms (85.56%) and identified 86,528 differentially expressed genes. We also identified 55 non-redundant full-length isoforms related to the flavonoid biosynthetic pathway. Pearson correlation analysis revealed that the expression levels of some key genes of the flavonoid biosynthesis pathway were significantly positively correlated with the flavonoid metabolites. Furthermore, we performed bioinformatics analysis (sequence and structural) of isoform_47029 (encoding flavanone 3-hydroxylase) and isoform_53692 (encoding flavonol synthase) to evaluate their potential biological functions. Finally, we validated gene expression levels of 12 flavonoid-related key enzyme genes using quantitative real-time PCR. Overall, this study provides full-length transcriptome information on G. longituba for the first time and valuable molecular resources for further research on the medicinal properties of this plant.

Comparative Analysis of Plastomes in Elsholtzieae: Phylogenetic Relationships and Potential Molecular Markers.

The Elsholtzieae, comprising ca. 7 genera and 70 species, is a small tribe of Lamiaceae (mint family). Members of Elsholtzieae are of high medicinal, aromatic, culinary, and ornamentals value. Despite the rich diversity and value of Elsholtzieae, few molecular markers or plastomes are available for phylogenetics. In the present study, we employed high-throughput sequencing to assemble two Mosla plastomes, M. dianthera and M. scabra, for the first time, and compared with other plastomes of Elsholtzieae. The plastomes of Elsholtzieae exhibited a quadripartite structure, ranging in size from 148,288 bp to 152,602 bp. Excepting the absence of the pseudogene rps19 in Elsholtzia densa, the exhaustive tally revealed the presence of 132 genes (113 unique genes). Among these, 85 protein-coding genes (CDS), 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (rps19 and ycf1) were annotated. Comparative analyses showed that the plastomes of these species have minor variations at the gene level. Notably, the E. eriostchya plastid genome exhibited increased GC content regions in the LSC and SSC, resulting in an increased overall GC content of the entire plastid genome. The E. densa plastid genome displayed modified boundaries due to inverted repeat (IR) contraction. The sequences of CDS and intergenic regions (IGS) with elevated variability were identified as potential molecular markers for taxonomic inquiries within Elsholtzieae. Phylogenetic analysis indicated that four genera formed monophyletic entities, with Mosla and Perilla forming a sister clade. This clade was, in turn, sister to Collinsonia, collectively forming a sister group to Elsholtzia. Both CDS, and CDS + IGS could construct a phylogenetic tree with stronger support. These findings facilitate species identification and DNA barcoding investigations in Elsholtzieae and provide a foundation for further exploration and resource utilization within this tribe.

Reassessing the genetic variability of Tectona grandis through high-throughput genotyping: Insights on its narrow genetic base.

Teak (Tectona grandis Linn. f.) is considered one of the most expensive hardwoods in the world. The dispersion of the species over the years has taken the teak beyond its first sources of diversity and little is known about the genetic origin and genetic variability. Thus, this study aimed to investigate the genetic diversity and genetic population structure existing in a representative teak germplasm bank collection. DNA was extracted from young leaves and each sample were genotyped by whole genome sequencing at 3 giga bases per sample, the sequences are aligned using the genome, and SNPcalls and quality control were made. To study the population structure of the genotypes, Bayesian variational inference was used via fastStructure, the phylogenetic tree was based on the modified Euclidean distance and the clustering by the UPGMA hierarchical method. Genetic diversity was analyzed based on the pairwise genetic divergence (Fst) of Weir and Cockerham. Genotyping by sequencing resulted in a database of approximately 1.4 million of variations SNPs were used for analysis. It was possible to identify four populations with considerable genetic variability between and within them. While the genetic variability in teak is generally known to be narrow, this study confirmed the presence of genetic variability scale in teak, which is contrary to what was initially expected.

The First Genome-Wide Mildew Locus O Genes Characterization in the Lamiaceae Plant Family.

Powdery mildew (PM) is a widespread plant disease that causes significant economic losses in thousands crops of temperate climates, including Lamiaceae species. Multiple scientific studies describe a peculiar form of PM-resistance associated at the inactivation of specific members of the Mildew Locus O (MLO) gene family, referred to as mlo-resistance. The characterization of Lamiaceae MLO genes, at the genomic level, would be a first step toward their potential use in breeding programs. We carried out a genome-wide characterization of the MLO gene family in 11 Lamiaceae species, providing a manual curated catalog of 324 MLO proteins. Evolutionary history and phylogenetic relationships were studied through maximum likelihood analysis and motif patter reconstruction. Our approach highlighted seven different clades diversified starting from an ancestral MLO domain pattern organized in 18 highly conserved motifs. In addition, 74 Lamiaceae putative PM susceptibility genes, clustering in clade V, were identified. Finally, we performed a codon-based evolutionary analysis, revealing a general high level of purifying selection in the eleven Lamiaceae MLO gene families, and the occurrence of few regions under diversifying selection in candidate susceptibility factors. The results of this work may help to address further biological questions concerning MLOs involved in PM susceptibility. In follow-up studies, it could be investigated whether the silencing or loss-of-function mutations in one or more of these candidate genes may lead to PM resistance.

Chloroplast Genome Structure and Phylogenetic Analysis of 13 Lamiaceae Plants in Tibet.

BACKGROUND: The chloroplast (cp) genome has unique and highly conserved characteristics and is therefore widely used in species identification and classification, as well as to improve the in-depth understanding of plant evolution. METHODS: In this study, the cp genomes of 13 Lamiaceae plants in the Tibet Autonomous Region of China were sequenced, assembled and annotated using bioinformatics methods. Phylogenetic trees were constructed to reveal the phylogenetic relationship of related species in the Lamiaceae. RESULTS: The results showed that all 13 cp genomes had a typical four-segment structure, including one large single-copy (LSC) region, one pair of inverted repeat (IR) regions and one small single-copy (SSC) region. The sequence lengths of the 13 cp genomes were between 149,081 bp and 152,312 bp, and the average GC content was 37.6%. These genomes contained 131-133 annotated genes, including 86-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. A total of 542 SSR loci were detected using MISA software. The repeat types were mostly single-nucleotide repeats, accounting for 61% of simple repeats. A total of 26,328-26,887 codons were detected in 13 cp genomes. According to the RSCU value analysis, the codons mostly ended with A/T. Analysis of IR boundaries showed that the other species were relatively conserved, except for Nepeta laevigata (D. Don) Hand.-Mazz., which differed in gene type and location on both sides of the boundary. By analysing nucleotide diversity, two highly mutated regions located in the LSC and SSC regions were identified in the 13 cp genomes. CONCLUSIONS: Using the cp genome of Lycium ruthenicum Murray as the outgroup, 97 cp genomes of the Lamiaceae were used to construct an Maximum Likehood (ML) phylogenetic tree, in which these species were divided into eight major clades, corresponding to eight subfamilies based on morphological classification. The phylogenetic results based on monophyletic relationships were consistent with the morphological classification status at the tribe level.

East Asian-North American disjunctions and phylogenetic relationships within subtribe Nepetinae (Lamiaceae).

Biogeographic disjunctions, including intercontinental disjunctions, are frequent across plant lineages and have been of considerable interest to biologists for centuries. Their study has been reinvigorated by molecular dating and associated comparative methods. One of the "classic" disjunction patterns is that between Eastern Asia and North America. It has been speculated that this pattern is the result of vicariance following the sundering of a widespread Acrto-Teritary flora. Subtribe Nepetinae in the mint family (Lamiaceae) is noteworthy because it contains three genera with this disjunction pattern: Agastache, Dracocephalum, and Meehania. These disjunctions are ostensibly the result of three separate events, allowing for concurrent testing of the tempo, origin, and type of each biogeographic event. Using four plastid and four nuclear markers, we estimated divergence times and analyzed the historical biogeography of Nepetinae, including comprehensive sampling of all major clades for the first time. We recover a well-supported and largely congruent phylogeny of Nepetinae between genomic compartments, although several cases of cyto-nuclear discordance are evident. We demonstrate that the three disjunctions are pseudo-congruent, with unidirectional movement from East Asia at slightly staggered times during the late Miocene and early Pliocene. With the possible exception of Meehania, we find that vicariance is likely the underlying driver of these disjunctions. The biogeographic history of Meehania in North America may be best explained by long-distance dispersal, but a more complete picture awaits deeper sampling of the nuclear genome and more advanced biogeographical models.

Mitochondrial Genome Sequence of Salvia officinalis (Lamiales: Lamiaceae) Suggests Diverse Genome Structures in Cogeneric Species and Finds the Stop Gain of Genes through RNA Editing Events.

Our previous study was the first to confirm that the predominant conformation of mitochondrial genome (mitogenome) sequence of Salvia species contains two circular chromosomes. To further understand the organization, variation, and evolution of Salvia mitogenomes, we characterized the mitogenome of Salvia officinalis. The mitogenome of S. officinalis was sequenced using Illumina short reads and Nanopore long reads and assembled using a hybrid assembly strategy. We found that the predominant conformation of the S. officinalis mitogenome also had two circular chromosomes that were 268,341 bp (MC1) and 39,827 bp (MC2) in length. The S. officinalis mitogenome encoded an angiosperm-typical set of 24 core genes, 9 variable genes, 3 rRNA genes, and 16 tRNA genes. We found many rearrangements of the Salvia mitogenome through inter- and intra-specific comparisons. A phylogenetic analysis of the coding sequences (CDs) of 26 common protein-coding genes (PCGs) of 11 Lamiales species and 2 outgroup taxa strongly indicated that the S. officinalis was a sister taxon to S. miltiorrhiza, consistent with the results obtained using concatenated CDs of common plastid genes. The mapping of RNA-seq data to the CDs of PCGs led to the identification of 451 C-to-U RNA editing sites from 31 PCGs of the S. officinalis mitogenome. Using PCR amplification and Sanger sequencing methods, we successfully validated 113 of the 126 RNA editing sites from 11 PCGs. The results of this study suggest that the predominant conformation of the S. officinalis mitogenome are two circular chromosomes, and the stop gain of rpl5 was found through RNA editing events of the Salvia mitogenome.

Evaluation of Candidate Reference Genes for Gene Expression Analysis in Wild Lamiophlomis rotata.

Lamiophlomis rotata (Benth.) Kudo is a perennial and unique medicinal plant of the Qinghai-Tibet Plateau. It has the effects of diminishing inflammation, activating blood circulation, removing blood stasis, reducing swelling, and relieving pain. However, thus far, reliable reference gene identifications have not been reported in wild L. rotata. In this study, we identified suitable reference genes for the analysis of gene expression related to the medicinal compound synthesis in wild L. rotata subjected to five different-altitude habitats. Based on the RNA-Seq data of wild L. rotata from five different regions, the stability of 15 candidate internal reference genes was analyzed using geNorm, NormFinder, BestKeeper, and RefFinder. TFIIS, EF-1α, and CYP22 were the most suitable internal reference genes in the leaves of L. rotata from different regions, while OBP, TFIIS, and CYP22 were the optimal reference genes in the roots of L. rotata. The reference genes identified here would be very useful for gene expression studies with different tissues in L. rotata from different habitats.

Conserved chloroplast genome sequences of the genus Clerodendrum Linn. (Lamiaceae) as a super-barcode.

BACKGROUND: The plants of the genus Clerodendrum L. have great potential for development as an ornamental and important herbal resource. There is no significant morphological difference among many species of the genus Clerodendrum, which will lead to confusion among the herbs of this genus and ultimately affect the quality of the herbs. The chloroplast genome will contribute to the development of new markers used for the identification and classification of species. METHODS AND RESULTS: Here, we obtained the complete chloroplast genome sequences of Clerodendrum chinense (Osbeck) Mabberley and Clerodendrum thomsoniae Balf.f. using the next generation DNA sequencing technology. The chloroplast genomes of the two species all encode a total of 112 unique genes, including 80 protein-coding, 28 tRNA, and four rRNA genes. A total of 44-42 simple sequence repeats, 19-16 tandem repeats and 44-44 scattered repetitive sequences were identified. Phylogenetic analyses showed that the nine Clerodendrum species were classified into two clades and together formed a monophyletic group. Selective pressure analyses of 77 protein-coding genes showed that there was no gene under positive selection in the Clerodendrum branch. Analyses of sequence divergence found two intergenic regions: trnH-GUG-psbA, nhdD-psaC, exhibiting a high degree of variations. Meanwhile, there was no hypervariable region identified in protein coding genes. However, the sequence identities of these two intergenic spacers (IGSs) are greater than 99% among some species, which will result in the two IGSs not being used to distinguish Clerodendrum species. Analysis of the structure at the LSC (Large single copy) /IR (Inverted repeat) and SSC (Small single copy)/IR boundary regions showed dynamic changes. The above results showed that the complete chloroplast genomes can be used as a super-barcode to identify these Clerodendrum species. The study lay the foundation for the understanding of the evolutionary process of the genus Clerodendrum.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory.

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

A highly contiguous genome assembly of red perilla (Perilla frutescens) domesticated in Japan.

Perilla frutescens (Lamiaceae) is an important herbal plant with hundreds of bioactive chemicals, among which perillaldehyde and rosmarinic acid are the two major bioactive compounds in the plant. The leaves of red perilla are used as traditional Kampo medicine or food ingredients. However, the medicinal and nutritional uses of this plant could be improved by enhancing the production of valuable metabolites through the manipulation of key enzymes or regulatory genes using genome editing technology. Here, we generated a high-quality genome assembly of red perilla domesticated in Japan. A near-complete chromosome-level assembly of P. frutescens was generated contigs with N50 of 41.5 Mb from PacBio HiFi reads. 99.2% of the assembly was anchored into 20 pseudochromosomes, among which seven pseudochromosomes consisted of one contig, while the rest consisted of less than six contigs. Gene annotation and prediction of the sequences successfully predicted 86,258 gene models, including 76,825 protein-coding genes. Further analysis showed that potential targets of genome editing for the engineering of anthocyanin pathways in P. frutescens are located on the late-stage pathways. Overall, our genome assembly could serve as a valuable reference for selecting target genes for genome editing of P. frutescens.

[Cloning of StHD1 and StHD8 from Schizonepeta tenuifolia and function of regulating glandular trichome development].

Hd-Zip, a unique transcription factor in plant kingdom, influences the growth, development, and secondary metabolism of plants. Hd-zip Ⅳ is thought to play an important role in trichome development of Schizonepeta tenuifolia. This study aims to explore the functions of StHD1 and StHD8 in Hd-zip Ⅳ subfamily in peltate glandular trichome development. To be specific, the expression patterns of the two genes and interaction between the proteins encoded by them were analyzed based on transcriptome sequencing and two-hybrid screening. The subcellular localization was performed and functions of the genes were verified in tobacco and S. tenuifolia. The results showed that StHD1 and StHD8 had high similarity to HD-Zip Ⅳ proteins of other plants and they all had the characteristic conserved domains of HD-Zip Ⅳ subfamily. They were located in the nucleus. The two genes mainly expressed in young tissues and spikes, and StHD1 and StHD8 proteins interacted with each other. The density and length of glandular trichomes increased significantly in tobacco plants with the overexpression of StHD1 and StHD8. Inhibiting the expression of StHD1 and StHD8 by VIGS(virus-induced gene silencing) in S. tenuifolia resulted in the reduction in the density of peltate glandular trichomes, the expression of key genes related to mono-terpene synthesis, and the relative content of limonene and pulegone, the main components of monoterpene. These results suggested that StHD1 and StHD8 of S. tenuifolia formed a complex to regulate glandular trichomes and affect the biosynthesis of monoterpenes.

Molecular and morphological survey of Lamiaceae species in converted landscapes in Sumatra.

Molecular biodiversity surveys have been increasingly applied in hyperdiverse tropical regions as an efficient tool for rapid species assessment of partially undiscovered fauna and flora. This is done by overcoming shortfalls in knowledge or availability of reproductive structures during the sampling period, which often represents a bottleneck for accurate specimens' identification. DNA sequencing technology is intensifying species discovery, and in combination with morphological identification, has been filling gaps in taxonomic knowledge and facilitating species inventories of tropical ecosystems. This study aimed to apply morphological taxonomy and DNA barcoding to assess the occurrence of Lamiaceae species in converted land-use systems (old-growth forest, jungle rubber, rubber, and oil palm) in Sumatra, Indonesia. In this species inventory, we detected 89 specimens of Lamiaceae from 18 species distributed in seven subfamilies from the Lamiaceae group. One third of the species identified in this study lacked sequences in the reference database for at least one of the markers used (matK, rbcL, and ITS). The three loci species-tree recovered a total of 12 out of the 18 species as monophyletic lineages and can be employed as a suitable approach for molecular species assignment in Lamiaceae. However, for taxa with a low level of interspecific genetic distance in the barcode regions used in this study, such as Vitex gamosepala Griff. and V. vestita Wall. ex Walp., or Callicarpa pentandra Roxb. and C. candidans (Burm.f.) Hochr., the use of traditional taxonomy remains indispensable. A change in species composition and decline in abundance is associated with an increase in land-use intensification at the family level (i.e., Lamiaceae), and this tendency might be constant across other plant families. For this reason, the maintenance of forest genetic resources needs to be considered for sustainable agricultural production, especially in hyperdiverse tropical regions. Additionally, with this change in species composition, accurate species identification throughout molecular assignments will become more important for conservation planning.

Mining, expression, and phylogenetic analysis of volatile terpenoid biosynthesis-related genes in different tissues of ten Elsholtzia species based on transcriptomic analysis.

We sequenced the leaf and inflorescence transcriptomes of 10 Elsholtzia species to mine genes related to the volatile terpenoid metabolic pathway. A total of 184.68 GB data and 1,231,162,678 clean reads were obtained from 20 Elsholtzia samples, and 333,848 unigenes with an average length of at least 1440 bp were obtained by Trinity assembly. KEGG pathway analysis showed that there were three pathways related to volatile terpene metabolism: terpenoid backbone biosynthesis (No. ko00900), monoterpenoid biosynthesis (No. ko00902), and sesquiterpenoid and triterpenoid biosynthesis (No. ko00909), with 437, 125, and 121 related unigenes, respectively. The essential oil content and composition in 20 Elsholtzia samples were determined by gas chromatography-mass spectrometry. The results showed that there were obvious interspecific differences among the 10 Elsholtzia species, but there were no significant differences between the different tissues among species. The expression levels of seven candidate genes involved in volatile terpenoid biosynthesis in Elsholtzia were further analyzed by quantitative real-time PCR. The results showed that HMGS had the highest expression among all genes, followed by GGPS4. In addition, there was not a significant correlation between the seven genes and the components with high essential oil contents. Combined with the essential oil components detected in this study, the possible biosynthetic pathway of the characteristic components in Elsholtzia plants was speculated to be a metabolic pathway with geraniol as the starting point and elsholtzione as the end product. Phylogenetic analysis was conducted using the nucleotide sequences of the geranyl diphosphate synthase candidate genes, and the results showed that genes related to the volatile terpenoid biosynthetic pathway may be more suitable gene fragments for resolving the Elsholtzia phylogeny.

Chromosome-level and haplotype-resolved genome provides insight into the tetraploid hybrid origin of patchouli.

Patchouli (Pogostemon cablin (Blanco) Benth.), a member of the Lamiaceae family, is an important aromatic plant that has been widely used in medicine and perfumery. Here, we report a 1.94 Gb chromosome-scale assembly of the patchouli genome (contig N50 = 7.97 Mb). The gene annotation reveals that tandem duplication of sesquiterpene biosynthetic genes may be a major contributor to the biosynthesis of patchouli bioactivity components. We further phase the genome into two distinct subgenomes (A and B), and identify a chromosome substitution event that have occurred between them. Further investigations show that a burst of universal LTR-RTs in the A subgenome lead to the divergence between two subgenomes. However, no significant subgenome dominance is detected. Finally, we track the evolutionary scenario of patchouli including whole genome tetraploidization, subgenome divergency, hybridization, and chromosome substitution, which are the key forces to determine the complexity of patchouli genome. Our work sheds light on the evolutionary history of patchouli and offers unprecedented genomic resources for fundamental patchouli research and elite germplasm development.

Sequence Analysis of the Complete Mitochondrial Genome of a Medicinal Plant, Vitex rotundifolia Linnaeus f. (Lamiales: Lamiaceae).

In the present study, we depicted the complete mitochondrial genome of a valuable medicinal plant, Vitex rotundifolia. The mitochondrial genome of V. rotundifolia, mapped as a circular molecule, spanned 380,980 bp in length and had a GC content of 45.54%. The complete genome contained 38 protein-coding genes, 19 transfer RNAs (tRNAs), and 3 ribosomal RNAs (rRNAs). We found that there were only 38.73% (147.54 kb), 36.28% (138.23 kb), and 52.22% (198.96 kb) of the homologous sequences in the mitochondrial genome of V. rotundifolia, as compared with the mitochondrial genomes of Scutellaria tsinyunensis, Boea hygrometrica, and Erythranthe lutea, respectively. A multipartite structure mediated by the homologous recombinations of the three direct repeats was found in the V. rotundifolia mitochondrial genome. The phylogenetic tree was built based on 10 species of Lamiales, using the maximum likelihood method. Moreover, this phylogenetic analysis is the first to present the evolutionary relationship of V. rotundifolia with the other species in Lamiales, based on the complete mitochondrial genome.

Whole Transcriptome Sequencing Unveils the Genomic Determinants of Putative Somaclonal Variation in Mint (Mentha L.).

Mint (Mentha L., Lamiaceae) is a strongly scented herb of the family Lamiaceae that is grown mostly by clonal propagation, making it a valuable species for the study of somaclonal variation and its phenotypic consequences. The recent introduction of a few species of mint in South America, followed by a presumably rampant propagation, make this region particularly ideal for studying the extent of somaclonal genetic diversity. Hence, the objective of this work was to offer a preliminary characterization of somaclonal genetically coding diversity of the mint in the northern Andes in order to address the question of whether somaclonal variants may have emerged despite relatively recent introductions in a region where mint is not native. A total of 29 clonally propagated specimens, collected in mint export farms in the province of Antioquia, a major region for mint production in the northwest Andes of Colombia, were genotyped using RNA sequencing (RNA-Seq). SNP calling was carried out from the leaves' transcriptome profiles of each plant by combining the GATK4 and TRINITY protocols, obtaining a total of 2033 loci across 912 transcripts with a minimum read depth of 20X and 4% of missing data. Unsupervised machine learning algorithms considered the K-means, AGNES and UPGMA approaches, all of which suggested three genetic clusters for M. spicata and a unique cluster for M. × piperita. The results indicate that at least two different origins of M. spicata reached the eastern region of the Antioquia province, clonally propagated in the locality ever since for local consumption and export. One of these ancestries had more population structure, possibly due to environmental or anthropological pressures that intervened in the fragmentation of this genetic group or to a higher somaclonal mutation rate. This work offers a first step into the study of the accumulation and transmission of presumably quasi-neutral somatic mutations at coding regions in an herbaceous clonally propagated scented species such as mint, likely favored by an expected population expansion after its Andean introduction. These ad hoc hypotheses warrant further study as part of future research.

Phylogeny-Aware Chemoinformatic Analysis of Chemical Diversity in Lamiaceae Enables Iridoid Pathway Assembly and Discovery of Aucubin Synthase.

Countless reports describe the isolation and structural characterization of natural products, yet this information remains disconnected and underutilized. Using a cheminformatics approach, we leverage the reported observations of iridoid glucosides with the known phylogeny of a large iridoid producing plant family (Lamiaceae) to generate a set of biosynthetic pathways that best explain the extant iridoid chemical diversity. We developed a pathway reconstruction algorithm that connects iridoid reports via reactions and prunes this solution space by considering phylogenetic relationships between genera. We formulate a model that emulates the evolution of iridoid glucosides to create a synthetic data set, used to select the parameters that would best reconstruct the pathways, and apply them to the iridoid data set to generate pathway hypotheses. These computationally generated pathways were then used as the basis by which to select and screen biosynthetic enzyme candidates. Our model was successfully applied to discover a cytochrome P450 enzyme from Callicarpa americana that catalyzes the oxidation of bartsioside to aucubin, predicted by our model despite neither molecule having been observed in the genus. We also demonstrate aucubin synthase activity in orthologues of Vitex agnus-castus, and the outgroup Paulownia tomentosa, further strengthening the hypothesis, enabled by our model, that the reaction was present in the ancestral biosynthetic pathway. This is the first systematic hypothesis on the epi-iridoid glucosides biosynthesis in 25 years and sets the stage for streamlined work on the iridoid pathway. This work highlights how curation and computational analysis of widely available structural data can facilitate hypothesis-based gene discovery.